Discrete Family of Dissipative Soliton Pairs in Mode-Locked Fiber Lasers 
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We numerically investigate the formation of soliton pairs (bound states) in mode-locked fiber 
ring lasers. In the distributed model (complex cubic-quintic Ginzburg-Landau equation) we observe 
a discrete family of soliton pairs with equidistantly increasing peak separation. This family was 
identified by two alternative numerical schemes and the bound state instability was disclosed by 
a linear stability analysis. Moreover, similar families of unstable bound state solutions have been 
found in a more realistic lumped laser model with an idealized saturable absorber (instantaneous 
response). We show that a stabilization of these bound states can be achieved when the finite 
relaxation time of the saturable absorber is taken into account. The domain of stability can be 
controlled by varying this relaxation time. 
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I. INTRODUCTION 

Since the experimental observation of stable solitary 
pulses (SPs) in optical fibers [l| their interactions at- 
tracted a great deal of interest. In the early stages the 
propagation of light pulses was described using the non- 
linear Schrodinger equation (NLSE). Relying on a per- 
turbation approach it was shown that SP interactions 
depend on the relative phase between them: in-phase 
solitons attract each other, while out-of-phase solitons 
repel and for an arbitrary phase the interac- 

tion appears more involved [H]- 

For a more realistic description of real systems dif- 
ferent perturbations were added to the NLSE, account- 
ing for losses [1] , intrapulse Raman scattering 0, [1, [l] , 
third-order dispersion 0, Q or spectral filtering. For 
the description of a complete communication line with 
a periodic arrangement of signal regeneration (semicon- 
ductor amplifier, saturable absorber) and propagation in 
the fiber (effects as nonlincarity and absorption) a dis- 
tributed description with the so-called complex cubic- 
quintic Ginzburg-Landau equation (CQGLE) was estab- 
lished PH . Another optical system of primary impor- 
tance, which is commonly modeled by this equation, is 
the mode-locked fiber ring laser, because it essentially 
consists of the same optical components in a closed loop 
PH . Despite its success in the description of mode- locked 
lasers, the assumptions used to derive the CQGLE are 
not always fulfilled in reality. 

In this contribution we investigate the formation of 
bound states in fiber lasers exploiting the CQGLE and a 
more realistic lumped model. From exact numerical sim- 
ulations of the CQGLE we obtain a discrete family of BS 
stationary solutions with different (equidistant) peak sep- 
aration. Moreover, we show that a lumped laser model 



* Electronic address: | aleksandr . zavyalov@uni- j ena.de] 



which accounts additionally for finite temporal effects of 
the absorber can have profound consequences for the BS 
stability. 

The CQGLE can be considered as a perturbed NLSE 
where the complex coefficients of dispersion and nonlin- 
carity account for spectral filtering and saturation, re- 
spectively, and a linear loss/gain is added. For certain 
parameter sets this equation was shown to support soli- 
tary wave solutions, called dissipative solitons, as well as 
bound states (BS) of pairs of SPs [H, [ll, [H, [lljp, [H, 
[Tgl . [20l . [2]| . Careful numerical investigations [la] have 
demonstrated BS formation from a pair of initially rest- 
ing SPs. Depending on the dispersion regime and the 
initial phase difference the BS profile can vary. Further- 
more it was shown [13, [III that out-of-phase and in-phase 
SPs can form unstable stationary BS solutions, because 
they represent saddle points in the phase plane. In the 
framework of the CQGLE stable BSs were discovered by 
Akhmediev et al. [2l| as stable two-soliton solutions with 
a 7r/2 phase difference of the peaks. They were observed 
for anomalous dispersion only. 

In an early work it was analytically shown that 
these dissipative equations exhibit a discrete family of 
two-soliton solutions with different (equidistant) peak 
separation. This analysis was based on both a perturba- 
tion approach for the NLS equation with small pumping 
and dissipative terms and the CQGLE with weak anoma- 
lous dispersion. To the best of our knowledge, however, 
such discrete families of BSs have not yet been found 
numerically or experimentally. 

In the first section of this paper we briefly present the 
equations used for modeling the fiber ring laser for dif- 
ferent approximations. We depart from the more gen- 
eral lumped description, where all individual elements in 
the cavity are modeled by separate equations and the 
pulse is sequentially propagated through these elements 
in the ring cavity. Then, we derive the CQGLE from 
the lumped model by averaging over the full cavity. We 
obtain definite relations between the coefficients of both 
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models which is essential for the comparison with exper- 
iments. 

The second section is devoted to the discrete family of 
stationary BS solutions derived from the CQGLE. In the 
numerical simulations we observe a discrete family of sta- 
tionary two-soliton solutions with different peak-to-peak 
separation which we identify as different BS levels. These 
results are obtained by means of two numerical schemes, 
viz., the solution of either the evolution or the stationary 
problem. Based on the results of the stationary analysis, 
a linear stability analysis is carried out where the insta- 
bility growth rates perfectly coincide with the results of 
the propagation model. 

For a more realistic description we take advantage of 
the lumped laser model in the third part, where the indi- 
vidual elements of the cavity are treated separately, and 
show that the different BS levels can be stabilized by ac- 
counting for a noninstantaneous, but very fast, response 
of the saturable absorber (SA), whereas for an instan- 
taneous absorber response the BSs remain unstable. In 
the latter case the SA is described by the commonly used 
instantaneous response approximation (ideal SA). 

Interestingly, by varying the SA response time and sat- 
uration level the BS solutions can be forced to change 
from stable to unstable behavior. Moreover, all these 
stable two-pulse solutions exhibit a multilevel nature and 
are located on a certain side of the phase plane depending 
on the level. 

In what follows we are restricting ourselves to the 
normal dispersion regime because this is the realm of 
the so-called all-normal-dispersion lasers, which attract 
a considerable deal of interest presently. Compared to 
other fiber lasers they have a simpler setup and allow 
for achieving higher pulse energies. Energies as high as 
256 nJ per pulse were achieved in a large-mode-area pho- 
tonic crystal fiber laser and more than 20 nJ for a 
usual single- mode fiber laser [2^. Pulses generated in 
all-normal-dispersion lasers depend nontrivially on the 
interplay of spectral filtering and self-phase-modulation 
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dH (1) 
i7|t/(z,t)p[/(z,i) , 



where U{z,t) is the envelope of the pulse, z is the prop- 
agation coordinate, t is the retarded time, /32 is the 
second-order dispersion (GVD) coefficient, and 7 repre- 
sents the fiber nonlinearity. g(z) is the saturable gain 
of the doped fiber and Ti is the dipole relaxation time 
(inverse linewidth of the parabolic gain). Assuming that 
the conditions are close to stationarity, the gain can be 
approximated by [131 
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where go is the small-signal gain, which is defined by the 
pumping level, and -E^j'" is the saturation energy. 

To describe the time-dependent semiconductor SA re- 
sponse we use the Agrawal/Olsson model [l^, which we 
term noninstantaneous SA response and is given by 
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where S{z,t) is the loss introduced by the absorber, 60 
is the small-signal loss, Treiax is the recovery time and 
Ef^ is the saturation energy. For Trciax < Tpuiso and 
taking into account that standard absorbers are thin, we 
obtain the well-known transmission equation for the SA 
in the instantaneous response approximation (ideal SA) 
from system (3) 
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II. MODELS 

For the numerical modeling we use a simple scheme 
of a ring fiber laser, which consists of a doped fiber, a 
saturable absorber and an output coupler. This laser 
model allows for including the dominant effects into the 
simulations and is still close to reality. In the lumped 
model the propagation through each element is treated 
separately. 

The propagation along the doped fiber is described 
by the modified nonlinear Schrodinger equation which 
is given by 0, [H, (when the carrier optical frequency 
equals the dopant's atomic resonance frequency) 



where Az is the length of the SA, SqAz defines the mod- 
ulation depth of the absorber and the saturation power 
is defined as P^at = Sf^/T^cia^. 

This lumped model describes the experimental setup 
appropriately but it is too involved for analytical stud- 
ies, including the linear stability analysis. In order to 
obtain a single equation with constant coefficients the 
periodic system is approximated by averaging over the 
full cavity length using the guiding-center soliton tech- 
nique [2^, [s^l, which relies on the field decomposition 
into a product containing a periodic part and an average 
amplitude. This model is justified if physically relevant 
changes appear upon several round-trips. For the sake of 
simplicity of the resulting equation the instantaneous re- 
sponse approximation for the SA (4) is used and the sat- 
uration of the fiber gain is neglected, i.e. we simplify Eq. 
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(2) to g{z) w go in the limit / \U{z,t)\^dt < Egf\ 

pulse 

Eventually we obtain a single evolution equation, which is 
a modified complex Ginzburg-Landau equation (MGLE) 
given by 
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where V is the normalized envelope of the pulse, Z = 
z/ Ld is the normalized averaged propagation coordinate, 
Ld = TQl\f5i\ is the dispersion length, r = t/To is the 
time normalized by the pulse duration Tq and the other 
parameters are defined as 
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1 - exp(-go£f) 
exp(ffo£f) - 1^ 

where Lf is the fiber length, which is assumed to be equal 
to the cavity length and k is the output loss. A further 
simplification can be achieved by a Taylor expansion of 
the last term of Eq. (4) up to second order in |V^P/^'s1vt''' ■ 
Ultimately we obtain the established CQGLE given by 



.dV Dd^V ,^,,2,, 

dZ 2 (7) 

= i5V + + ie\V\^V - V , 

where the linear and nonlinear gain/loss coefficients arc 
defined as 



p, e 
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Hence, we can unambiguously relate the coefficients of 
the CQGLE and MGLE to the coefficients of the lumped 
model which reflect the experimental laser parameters. 



III. THE CQGLE MODEL 

A. Bound state solutions 

To get a general picture of the two-SP interactions in 
the fiber laser we start our considerations from the dis- 
tributed laser model, CQGLE (7). The reason is twofold. 



first most previous papers rely on this model and second, 
it is more appropriate for identifying the stationary so- 
lutions and applying the linear stability analysis. 

From previous investigations [H, [li, [H, [ll, [13, [H, 
[TqI . [20| . [2i| it is known that BS formation is caused by 
linear and nonlinear dissipative effects (right hand side 
of Eq. (7)). Thus to simplify the analysis the number of 
system parameters was reduced from four (5, /3, e, fi) to 
two {j3, k) such as [l6| : 



5 = -fc/3/3, e = (}/2 + k(i, p = (5/8) fc/3, (9) 

where the parameter k characterizes the deviation from 
the particular point (fc = 0) where arbitrary-amplitude 
solutions exist |3ll |. 

We are working in the normal dispersion regime D = 
— 1 with relatively strong spectral filtering (3 = 1. To be 
sufficiently far from the singularity we choose k = 0.5 
where stable SPs exist. As initial conditions wc use an 
in-phasc superposition of two resting pulses separated by 
a certain distance. In the framework of the CGQLE the 
initial pulse profile is a critical issue with respect to BS 
formation. Thus for a better convergence we use as input 
two SPs, which are numerically obtained stationary one- 
soliton solutions for the respective parameter set . 

To get a more complete picture of two-pulse solutions 
we consider the evolution for a wider range of initial sepa- 
rations and longer propagation distance than in previous 
works, see Fig. 1. We observe stationary BS formation 
beyond a certain distance (similar to [171) ^^^^ its suc- 
ceeding destabilization, i.e. fusion of the pulses (below 
we discuss this instability in more detail). For an initial 
separation exceeding a critical value the SPs repel each 
other, was also observed in [H, |T3]. However, a more 
careful study discloses that this separation is no ordi- 
nary repulsion but the formation of another BS with a 
larger peak separation (see Fig. 1). In previous works 
it was usually assumed that the interaction becomes too 
weak for such large peak separations and, hence, this sce- 
nario was not thoroughly investigated. The evolution to 
the second BS level has a much slower dynamics than 
that for the first one and needs an approximately 50 
times longer propagation distance to form the station- 
ary solution. This difference arises from the weaker SP 
interaction. Except for the larger formation distance and 
the slower formation dynamics there is no principal dif- 
ference between these two BS levels. Similarly, beyond 
another critical initial separation we have observed the 
evolution towards a third BS level. Unfortunately, these 
calculations are too time-consuming, but intuitively it is 
clear that even higher order BS levels should occur. Evi- 
dently, it may be anticipated that the dynamics becomes 
progressively slower for higher levels and that all of them 
are unstable. 

To double-check these expectations and to identify ar- 
bitrary BS levels, we solve the stationary problem using 
the relaxation and Newton's methods in combination. 
Figure 1 (right) shows results obtained from both the 
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FIG. 1: Evolution of the separation for initially in-phase 
pulses with different initial separation. The levels of station- 
ary in-phase BS solutions are displayed on the right axis on 
the right side. The peak separation is given in units of r; 
model: CQGLE (7); parameters: k — 0.5, /3 = 1. 
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FIG. 2: (Color online) Intensity and phase profiles for the first 
three levels of the (a) in-phase and (b) out-of-phase stationary 
BS solutions; model: CQGLE (7); parameters: k — 0.5, /? = 
1. 



propagation and from the stationary analysis, which are 
in perfect agreement in the intermediate section where 
the propagation is stationary. The solution of the sta- 
tionary problem yielded six levels of in-phase BSs with 
different peak separation. It is interesting to see that 
the peak separation S changes from level to level by 
a constant value (approximately AS' = 6 for the cur- 
rent set of parameters) and obeys the simple equation 
'^in-phaso ~ 2.5 + 6n whcrc n — 1, 2, 3, ... designates the 
BS level. 

Also the first six levels of the out-of-phase stationary 
BS solutions were found. Similarly, the dependence of 
the peak separation on the level (except the first one) 
can be approximated by S"^^_^f_p^^^ w 5.5 -I- 6n, where 
n ~ 2, 3, 4, while for the first level Sl^^_^f_py^^g^ = 4.7. 
We guess that this peculiarity of the first level is caused 
by the strong SP overlap resulting in enhanced nonlinear 
effects. 

We note that for a 7r/2 phase difference between the 
pulses neither stable nor unstable BS solutions have been 
found in the normal dispersion regime. 

To understand the multilevel nature we consider the 
BS intensity and phase profiles for the first three levels, 
shown in Fig. 2. From this figure we may draw two con- 
clusions. First, the intensity profile of in-phase (out-of- 
phase) BS looks always the same, independently of the 
level, and has a characteristic small dip in the center. 
Second, there is a definite phase relation between the BS 
constituents and between the levels. For example, for in- 
phase BS levels the phase difference between the phase 
maximum and the phase minimum (at the center) is a 
multiple of tt and can be written as A(/j?Jj_pjjg^gg sa nn, 
where n = 1,2,3,... is the BS level. For out-of-phase 
solutions this phase difference can be approximated by 
^(^out- of -phase ~ 0-8 + 7rn, where n = 2,3,4,..., while 
for the first level A^out-of-phase = 3-6. 

Physically, the equidistant multilevel nature of the BSs 
can be understood as constructive interference between 
SPs and can be clearly seen from the almost linearly 



growing phase of the tails (Fig. 2). On the other side, 
our results are in qualitative agreement with Malomed's 
analysis [ij, [l^ where the oscillating tails of the SPs 
(linear phase change) evoke a periodic interaction poten- 
tial against the pulse separation, in spite of the fact that 
this analysis [ij, [l^ is based on a NLSE perturbation 
approach and is valid just for the anomalous dispersion 
regime. Moreover, we propose the magnitude Atp as a 
convenient parameter for the definition of the BS level in 
simulations or experiments. 

We visualize the evolution of both pulses in a phase 
plane (p{z), 4){z)) where p{z) is the peak separation and 
(j>{z) their relative phase (difference of the phase between 
peaks), see Fig. 3. In this plot we analyze the evolution 
trajectories of two pulses for two relevant cases. In the 
first case we consider two out-of-phase SPs with an initial 
separation of 5 = 7 which are evolving towards the first 
level out-of-phase stationary BS solution. In the second 
case the evolution of a pair of in-phase SPs with a larger 
initial separation oi S = 16.1 (evolving towards the sec- 
ond level in-phase stationary BS solution) is displayed. 
The evolution trajectory for the first case, which may be 
partly described by a circle indicated by '1' in Fig. 3, 
is well studied The radius of this circle equals the 
distance between the peaks of the first level out-of-phase 
BS. 

To date, it was not shown that in the second case, 
where the initial separation between the solitons is 16.1, 
the evolution trajectory partly consists of four circles. 
From the plot (Fig. 3) we can recognize that each circle 
corresponds to a certain BS. The radius of the biggest 
one equals the distance between pulses of the second level 
in-phase BS, being in agreement with the previous case. 
The other circles '3', '2' and '1' correspond to the second 
level out-of-phase, first level in-phase and first level out- 
of-phase stationary BS solutions, respectively. 

It is interesting to note that the evolution trajectory 
of two in-phase or out-of-phase SPs attains simple geo- 
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FIG. 3: (Color online) Soliton trajectories in the phase plane 
for: two initially out-of-phase solitons with separation 7 (line 
A) and two initially in- phase solitons with separation 16.1 
(line B). IL - first level stationary in-phase (red right) and 
out-of-phase (blue left) bound state solution, 2L - second level 
stationary in-phase (red right) and out-of-phase (blue left) 
bound state solution; model: CQGLE (7); parameters: k = 
0.5, 13 = 1. 



metrical forms. Moreover, during the evolution they pass 
through all possible in-phase and out-of-phase stationary 
BS solutions with ever smaller peak separation. 

This can be explained by the earlier finding that the 
phase difference between the pulses depends on their sep- 
aration as a consequence of interference phenomena be- 
tween them. 



FIG. 4: (Color online) Bound state instability. Two in- 
phase solitons form a quasi-stationary bound state but fuse 
asymmetrically eventually; model: CQGLE (7); parameters: 
k = 0.5, (3 = 1. 
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B. Linear stability analysis 

In previous papers [13, [H, HH it was shown that in- 
phase and out-of-phase first level BS are unstable. From 
our propagation simulations we anticipate that the corre- 
sponding higher level BSs arc unstable as well. In order 
to confirm this hypothesis and to provide a more com- 
plete picture we carry out a linear stability analysis of the 
higher BS levels. We start with the first level in-phase 
stationary solution. Its formation, the typical instability 
behavior and final fusion of the two-pulse excitation are 
shown in Fig. 4. It is interesting to note that the insta- 
bility exhibits a temporal asymmetry in the absolutely 
symmetric CQGLE. 

To perform the linear stability analysis we 
perturb the stationary solution Vo(i) exp(ifcoz) 
as V{z,t) = Vo{t) cxp{ikoz) + 5V{z,t) = 
[Vo{t) + f exp{Xz) + geyip{X*z)]exp{ikoz), where Vo{t) 
was numerically calculated and /, g are the small 
amplitudes of the perturbation. The substitution of 
this ansatz into the CQGLE (7) and its respective 
linearization in /, g leads to a system of linear, homo- 
geneous equations. The solvability condition provides 
the eigenvalues A. Eigenvalues with positive real part 



FIG. 5: (Color online) Complex eigenvalues for both the first 
level unstable stationary in-phase BS and stable SP solu- 
tion. The arrows denote the directions where the eigenvalues 
are shifted for higher levels BSs with decreasing interaction; 
model: CQGLE (7); parameters: k = 0.5, (3 = 1. 



indicate that the perturbation grows upon propagation 
and the stationary solution destabilizes. In Fig. 5 the 
calculated eigenvalues for the unstable first level in-phase 
BS are displayed and compared with the eigenvalues of 
the stable SP solution. 

Evidently, in the limit of two infinitely separated SPs 
the eigenvalues coincide with those of the stable single 
SP solution. In this case the zero eigenvalue exhibits a 
fourth order degeneracy corresponding to the two phase 
and two translation modes. If the peak separation de- 
creases the pulses start to perturb each other affecting 
the eigenvalues. Because of interaction two eigenvalues 
branch off from the zero eigenvalue and move in oppo- 
site directions on the real axis for decreasing separation. 
Thus, the eigenvalue with the positive real part is respon- 
sible for the instability of the BS solution, see Fig. 5. The 
two eigenmodes corresponding to these eigenvalues are 
shown in Fig. 6. They are formed as linear combinations 
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TABLE I: Re(A) for different levels of stationary BS solutions; 
model: CQGLE (7); parameters: k = 0.5, /3 = 1. 
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FIG. 6: (Color online) Left: Perturbation modes (dotted blue 
line and dashed red line) of the first level in-phase BS (solid 
black line). Right: Evolution trajectories of points A and B 
of the stationary BS solution. The amplitude of the pertur- 
bation modes is given in arbitrary units; model: CQGLE (7); 
parameters: k = 0.5, /3 = 1. 



of the four individual eigenmodes of the original SP (two 
phase modes, two translation modes) shifted by the dis- 
tance S. The peaks of both modes arc shifted towards 
the center and do not coincide with the peaks of the BS. 
Also because of SPs interaction another four eigenvalues 
are symmetrically split off the spectral continuum and 
move away from it for decreasing BS level. However, all 
of them have a negative real part and do not affect the 
stability of the bound state solution. 

From the results presented in Figs. 5 and 6 we may 
conclude that the antisymmetric mode is responsible for 
the BS instability. In the propagation simulations this 
mode is excited by noise and then grows with Re(A) > 0. 
To confirm the results of the linear stability analysis we 
add white noise to the stationary BS solution and plot 
the evolution trajectories of two points (fixed on the time 
coordinate, calculated in the propagation scheme). The 
points are chosen to be close to the amplitude maxima of 
the unstable mode, see Fig. 6 (left). The growing anti- 
symmetric eigenmode evokes an increase of the amplitude 
in one point and a decay in the other point upon prop- 
agation. Figure 6 (right) shows that the results of the 
linear stability analysis and the exact numerical propa- 
gation simulations perfectly coincide up to the distance 
where the linear stability analysis ceases to be valid. 

As expected, Re(A) calculated for other BS levels is 
smaller for higher levels, sec Table HI This dependency 
correlates well with the slower dynamics (evolution to the 
BS solution and instability dynamics) discussed in Fig. 
1. From a mathematical and practical point of view the 
growth rate Re(A) for higher BS levels can be very close 
to zero, but physically it remains an unstable BS. 

We varied the CQGLE parameters /3, k and 5 to find 
in-phase or out-of-phase stable BSs but failed. 
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IV. THE LUMPED MODEL 

In this section we study BS formation in fiber ring 
lasers using the more realistic lumped model and exper- 
imentally accessible laser parameters. The aim of this 
study is manifold. First of all, this model reflects bet- 
ter the experimental situation. So, it is interesting to 
compare the results with those of the distributed model 
which is easier to handle but contains severe approxima- 
tions. Moreover, we will study the effect of the finite 
SA response time (noninstantaneous model) on the BS 
stability. A stabilization has been proven in the spatial 
domain by taking advantage of nonlocality (see |32l | and 
references therein). Here, the noninstantaneous response 
plays the role of nonlocality in the spatial domain. Even- 
tually we will summarize the differences between the re- 
sults provided by the lumped and the distributed model 
in order to derive the limits of validity of the latter. 

For the simulation we use parameters closely related to 
previous experiments at a carrier wavelength of 1030 nm 
[13, m, [13 ■ The length of the absorber and the small- 
signal loss were adjusted to a modulation depth of 30%, 
the relaxation time of the fast, but noninstantaneous, SA 
is 500 fs, its saturation energy amounts to 16.7 pJ and 
the respective saturation power to Pgat = 33.4 W. The 
output loss is equal to 30%. For the doped fiber we as- 
sume if = 1 m, /32 = 0.024 ps^m-i, 7 0.005 W-^m.-^ 
and i?^t™ = 1 nJ. Amplitude and phase profiles were 
always recorded after the output coupler. 

In order to trigger BS formation we varied the filter 
strength, corresponding to the inverse gain bandwidth 
Ti of the laser material, and the small-signal gain. 



A. Instantaneous SA approximation 

In a first step we aim at reproducing the BS levels ob- 
tained in the CQGLE model. For this purpose we start 
with the instantaneous SA approximation, Eq. (4). As 
initial condition we use two in-phase or out-of-phase rest- 
ing small amplitude Gaussian pulses which converge to 
the BS. Because of the additional stabilization effect of 
the saturated gain the current model is less critical with 
respect to the initial pulse profiles than the CQGLE. We 
observe BS formation for Ti = 300 fs and go = 0.75 m~^. 
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FIG. 7: (Color online) Intensity and phase profiles for first 
(solid line 1) and second (dashed line 2) levels of the stationary 
BS solutions; (a) - in-phase, (b) - out-of-phase; lumped model 
(parameters, see text), instantaneous SA approximation (4). 
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FIG. 9: (Color online) Pulse trajectories in the phase plane for 
two pulses (a) initially in-phase with separation 6.38 ps and 
(b) initially out-of-phase with separation 5.25 ps. IL - first 
level BS, 2L - second level BS; lumped model (parameters, 
see text), instantaneous SA approximation (4). 



To visualize the interactions between both pulses we 
again take advantage of the phase plane. We consider 
the evolution trajectories for initially in-phase and out- 
of-phase pulses, see Fig. 9. For in-phase pulses the tra- 
jectory is partly described by four circles. As we already 
know, each circle (i?i, i?2, -R3, ^4) corresponds to a 
certain level of out-of-phase or in-phase stationary BS 
solutions, shown in Fig. 9. There is essentially no differ- 
ence between the evolution trajectories provided by the 
lumped and the distributed model. 
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FIG. 8: (Color online) Bound state instability. The two in- 
phase pulses form the stationary bound state and fuse asym- 
metrically eventually; lumped model (parameters, see text), 
instantaneous SA approximation (4). 



The first and second stationary in-phase BS levels are 
displayed in Fig. 7. In-phase BSs exhibit a small charac- 
teristic central peak in the intensity profile and the dif- 
ference between phase maximum and minimum (at the 
center) is a multiple of ~ tt. In agreement with our previ- 
ous results the formation of the second level in-phase BS 
needs approximately 80 times more (40,000) round trips 
than the 1st level (500). However, also these stationary 
BS solutions are unstable with an asymmetric instability 
behavior which is similar to the CQGLE BSs, see Fig. 8. 

For out-of-phase initial conditions we have also ob- 
served first and second level stationary out-of-phase BSs, 
see Fig. 7. They exhibit a similar formation dynamics 
and instability behavior as in-phase solutions. In par- 
ticular, the formation of a level 2 out-of-phase BS needs 
approximately 30 times more (3,000) round-trips than 
that for a level 1 out-of-phase BS (100). 



B. Differences between models 

For the sake of comparison between the different mod- 
els and of evaluating their validity we use the instan- 
taneous SA approximation and compare results provided 
by the lumped model, which serves as a benchmark, with 
both distributed models, the MGLE (5), which accounts 
for the full SA saturation, and the CQGLE (7), which 
relies on the Taylor expansion of the SA saturation. 

Figure 10 shows stable SP solutions obtained by both 
the lumped and the MGLE model. For the lumped mod- 
eling the saturation energy of the gain was reduced by 
a factor of two, because we consider just a single pulse. 
The MGLE parameters were calculated from those of the 
lumped model, except that for the distributed model a 
proper magnitude of the small-signal gain (6) was ad- 
justed as go = 0.6388 m^^, corresponding to the gain in 
the lumped model for stationary conditions. It is inter- 
esting to note that for the relevant parameter set stable 
SP solutions exist in the MGLE, but not in the CQGLE 
model. This is a clear indication that the saturation be- 
havior plays a pivotal role. Usually mode-locked lasers 
are working in a regime where the absorber is saturated, 
i.e., iFlV-PsaT" < 1 is not satisfied and the CQGLE is 
incorrect. Hence, one must not expect solutions for simi- 
lar parameters in the CQGLE model on the one side and 
the lumped (or MGLE) model on the other side. Never- 
theless, the MGLE and the lumped model with an ideal 
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FIG. 10: (Color online) Comparison between the lumped and 
the MGLE model at the same parameters set. (a) - intensity 
profile of stable SP solutions, (b) - evolution trajectories of 
two pulses in the phase plane; instantaneous SA approxima- 
tion (4). 



SA show quite a good agreement regarding the stable SP 
solutions, see Fig. 10. 

Further on we consider the evolution of two SPs in the 
phase plane, see Fig. 10 (in this case wc use the original 
value of the saturation energy fi'^t'" = 1 nJ). Although 
minor differences occur the evolution trajectories of the 
two pulses behave quite similar for both the MGLE and 
the lumped model. 



C. Noninstantaneous SA response 



6.0 

5.5 
I 5.0 
I 4.5 

CD 
CD 

g-4.0 
w 

I 3-5 

CD 
O 

■^3.0 
2.5 



(a) 



-T ,= 

relax 

T , = 150 fs- 



I , - OOU IS. 



V i/ V V 




(b) 



---T 

re 



= 

= 150fs 
: 350 fs 



6.0 
5.5 
5.0 
4.5 
4.0 
3.5 
3.0 
2.5 







3000 



6000 



9000 



3000 



6000 



9000 



round-trip number 



round-trip number 



FIG. 11: (Color online) Bound state stabilization. Distance 
between two pulses as a function of round-trips for initially 
(a) - in-phase and (b) - out-of-phase pulses for different SA 
relaxation times. In the limit Trciax = 0, pulses build (a) the 
level 1 in-phase stationary BS solution, (b) the level 2 out-of- 
phase stationary BS solution; lumped model (parameters, see 
text), noninstantaneous SA response (3). 




In a next step, relying on the lumped model, we are 
attempting to achieve BS stabilization by accounting for 
the noninstantaneous SA response (3). This model re- 
flects the experimental situation because typical SA re- 
sponse times vary between 300 fs and 12 ps, which fre- 
quently compares to the pulse length. We use the pa- 
rameters from the previous section A, in order to keep 
the saturation power constant, and account for the vari- 
ation of the relaxation time by a proper change of the 
saturation energy Ef^. 

In these simulations as initial condition we use two 
resting pulses with small amplitudes. At the beginning 
the pulses propagate 100 round trips in the ideal SA ap- 
proximation and just then we switch on the noninstan- 
taneous response. This procedure is required to arrive 
at a definite initial condition for large amplitude bound 
states. 

Trivially, for a vanishing relaxation time the system ex- 
hibits unstable stationary BS solutions identical to them 
of the instantaneous SA approximation. With increasing 
relaxation time (about 150 fs) we observe a tendency to- 
wards BS stabilization with moderate oscillations. Even- 
tually the BSs stabilize for a relaxation time (about 350 
fs) close to the experimental situation, see Figs. 11, 12. 
With a further increase of the relaxation time the arising 
stable BSs exhibit a fixed peak separation. The solu- 
tions move backwards in the reference frame due to the 
temporal effects in the absorber, see Fig. 12. 




500 1000 1500 2000 2500 3000 3500 4000 

round-trip number 

FIG. 12: (Color online) Oscillating BS solution for Trdax = 
150 fs; lumped model (parameters, see text), noninstanta- 
neous SA response (3). 



It is interesting to note that the resulting stationary 
and oscillating stable BS solutions do not depend on the 
initial phase difference between the pulses, see Fig. 11. 
This is in stark contrast to the CQGLE model, where 
the initial phase difference defines the nature of the sta- 
tionary BS solution. There, two in-phase or out-of-phase 
pulses with identical separation evolve towards different 
levels (see the curve for Trdax = in Fig. 11). More- 
over, we have observed stationary and oscillating stable 
BS solutions of the next level with a larger peak separa- 
tion. To obtain them the initial distance between pulses 
was increased and proper figures for the relaxation time 
were chosen. Higher level BSs need approximately 100 
times more round-trips for their formation, which is in 
agreement with previous results. 

The nature of stationary and oscillating stable BS so- 
lutions can be understood by inspecting the trajectories 
in the phase plane, see Fig. 13. We evaluate the two 
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FIG. 13: (Color online) Trajectories showing the two-pulse 
evolution; a) 1st level stable BS solution Troiax = 350 fs (blue 
line A), 2nd level stable BS solution T^oiax = 75 fs (red line B) 
and (b): 1st level oscillating BS solution Trdax = 150 fs (blue 
line A), 2nd level oscillating BS solution Treiax = 5 fs (red line 
B); lumped model (parameters, see text), noninstantaneous 
SA response (3). 



pulse dynamics in phase space for four different cases. 
To generate the first and second level BSs we use two 
pulses with an initial separation of 4.9 ps and 6 ps, re- 
spectively. Moreover, for each case we adjust the proper 
relaxation time in order to obtain a stationary or an os- 
cillating stable BS solution. Figure 13 shows that sta- 
tionary BS solutions are fixed points in the phase plane 
which are located near = 7r/2. Oscillating solutions 
move either in the lower or upper half plane between two 
adjacent circles. 

The conclusion to be drawn is that stationary stable BS 
solutions (Fig. 14) of the lumped model with noninstan- 
taneous SA are neither in- nor out-of-phasc but exhibit a 
phase difference of about tt/2 between the pulses. Imme- 
diately the question arises whether there are stable BS 
solutions with tt/2 phase difference between the pulses 
in the instantaneous models (CQGLE or lumped model 
with instantaneous SA). In the lumped model no stable 
or even stationary BS solutions with that phase differ- 
ence could be identified. For the CQGLE model, but 
just for anomalous dispersion, this case has been studied 
earlier [2l|, and it was shown, that stable BS solutions 
can exist. 

Nevertheless, stable BS solutions in the lumped model 
with noninstananeous SA differ considerably from those 
discussed in [2l|. First, they exist for a completely dif- 
ferent set of parameters (even opposite dispersion) and 
second, their behavior in phase space (Fig . 13) exhibits 



a large difference to those studied in [2lJ. Because all 
SAs have a finite relaxation time the conclusion can be 
drawn that the parameter range for the existence of BSs 
is much wider than assumed in [2l| . 

As already mentioned BS stabilization is evoked by the 
noninstantaneous SA response in analogy to the spatial 
case where a nonlocality stabilizes localized solutions [1^ . 
Here, the retarded temporal SA response is asymmetric 
in time, leading to the minor deviation of the fixed point 
from (j) = 7r/2. By contrast, the stable BS solutions of 






time (ps) 

FIG. 14: (Color online) The intensity and phase profiles for 
1st and 2nd level stable BS solutions calculated for Trciax = 
350 fs (blue solid line 1), T^oiax ~ 75 fs (red dashed line 2). 
The phase differences are Ayi = 5.1 and A(^2 ~ 6.6; lumped 
model (parameters, see text), noninstantaneous SA response 
(3). 
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FIG. 15: Existence domains for stable (stationary and os- 
cillating) bound state solutions; (a) first level solution, (b) 
second level solution; lumped model (parameters, see text), 
noninstantaneous SA response (3). 



the CQGLE appear exactly for = 7r/2 |21| . 

The region of existence of stable BSs is displayed 
in Fig. 15. It can be recognized that the regions for 
two neighboring BS levels differ, especially the bottom 
boundaries. This means that oscillating solutions of the 
first and second level can be observed for different relax- 
ation times and the same small-signal gain. The open 
boundary indicates that the bottom margin is very close 
to the zero axis. 

A distinctive feature of stable BS solutions is their 
behavior in the phase plane. The two initially resting 
pulses, for a certain range of initial separation, always 
evolve towards a single fixed point independently of the 
initial phase difference. This is in stark contrast to the 
stable BS solutions studied in Ref. [2l| where two sim- 
ilar solutions are possible, lying symmetrically on oppo- 
site sides of the phase plane. Moreover, the neighboring 
levels of BSs arc situated in different half spaces, see Fig. 
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13. This is due to the lack of invariancc of the time scale 
if relaxation matters. Thus leading and trailing pulses 
arc not equivalent. 

Figure 13 shows that each fixed point on the phase 
plane lies between two neighboring circles. Thus one 
might anticipate that there is a fixed point between circle 
1 and 2 too. But in our simulations we did not succeed 
in finding any stable BS in that domain. 

V. CONCLUSION 

In conclusion, we have theoretically and numerically 
investigated the formation and stability of BSs in a mode- 
locked fiber ring laser with SA in the normal disper- 
sion regime. We started our investigations with the dis- 
tributed description of the laser which is given by the 
commonly used CQGLE. For the first time, to the best 
of our knowledge, we have observed a discrete family of 
stationary BS solutions with different peak separation. 
Moreover, the observed BS levels exhibit an equidistantly 
increasing peak separation with a well-defined phase re- 
lation. All levels of in-phase and out-of-phase BSs turned 
out to be unstable for the investigated system parame- 
ters. A linear stability analysis has shown that this insta- 
bility is caused by the antisymmetric perturbation mode. 
We found that both the formation dynamics of BSs and 
their succeeding decay is slower for higher levels because 
of the weaker interaction. The evolution trajectories in 
the phase plane of two SPs with relatively large initial 
separation consist of simple geometrical forms. 

All these findings could be confirmed by using a more 



realistic lumped laser model but maintaining the instan- 
taneous SA response. 

Stabilization of the BSs was achieved in the lumped 
laser model by taking the fast, nevertheless noninstanta- 
neous SA response into account. Using the SA relaxation 
time as control parameter a continuous transition from 
unstable to stable BS solutions may be achieved passing 
a domain of oscillating BS solutions. These stationary 
and oscillating stable BSs have a discrete multilevel na- 
ture too. The distinctive feature of these BS solutions 
is expressed in the phase plane. A stable BS solution 
of a certain level represents a fixed point in the phase 
plane, which lies in the prefered half space. For oscillat- 
ing BSs the evolution is described by circulations between 
two neighboring circles in either the lower or upper half 
space depending on the level. The stabilization mecha- 
nism identified shows that all lasers with semiconductor 
absorber have a wider region of stable BS solutions than 
previously assumed [2l] |. Moreover, we have shown that 
a discrete family of stable or oscillating BS solutions with 
equidistant peak separation can be obtained in such laser 
systems. 

The results presented can be used in optics commu- 
nication lines for generating multipulse trains with well 
defined separation distances or in fiber lasers for an ir- 
regular or regular pulse array generation. 
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